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We analyze correlations between the event-plane angles in different intervals of pseudorapidity 
within the 3+1-dimensional viscous hydrodynamics with the Glauber-model initial conditions. As 
predicted earlier, the fluctuations in the particle production mechanism in the earliest stage, together 
with asymmetry of the emission profiles in pseudorapidity from the forward- and backward going 
wounded nucleons, lead to the torque effect, namely, decorrelation of the event-plane angles in 
distant pseudorapidity bins. We use two- or three-bin measures of correlation functions to quantify 
the effect, with the latter compared to the recent data from the CMS collaboration. We find a 
sizable torque effect, with magnitude larger at RHIC than at the LHC. 

PACS numbers: 25.75.-q, 25.75Gz, 25.75.Ld 


I. INTRODUCTION 

Fluctuations of the shape of the fireball created in rel¬ 
ativistic heavy-ion collisions give rise to event-by-event 
fluctuations of the collective flow. The main examples 
of such phenomena are the elliptic flow fluctuations [T], 
the triangular flow |2j, the mean transverse momentum 
fluctuations [3], the broad event-by-event distribution of 
the flow coefficients [4], or breaking of the factorization 
of the two-particle correlations [5j. The fluctuations of 
the initial distribution of sources give rise, due to the hy¬ 
drodynamic expansion, to event-by-event fluctuations of 
the collective flow }6]. 

A related aspect is the rapidity dependence of the col¬ 
lective flow. Two-particle correlations in relative pseu¬ 
dorapidity and relative azimuthal angle exhibit ridge-like 
structures extended in the longitudinal direction. These 
long-range fluctuations in rapidity indicate the presence 
of longitudinal correlations in the initial stage of the 
collision. The approximate preservation of the trans¬ 
verse shape of the fireball along the longitudinal direction 
yields strong correlations of the flow event-plane angles 
at different rapidities, which allows one to extract the 
flow coefficients using the event plane direction defined 
in a forward rapidity interval. 

A mechanism leading to decorrelation of event-planes 
at different rapidities, the torque effect, has been pro¬ 
posed in Ref. [7j. One expects that at forward (backward) 
rapidities the event-plane direction is predominantly de¬ 
termined by the participant nucleons from the forward 
(backward) going wounded nucleus. While strongly cor¬ 
related, the event-plane directions defined by the forward 
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(backward) going wounded nucleons are different due to 
event-by-event fluctuations. At intermediate rapidities 
the shape of the interaction region should interpolate be¬ 
tween the forward and backward limits. As a result, the 
orientation of the event-plane changes along pseudorapid¬ 
ity in a given event: the torque bends to the left or right, 
with a fluctuating angle. The details and the strength 
of the torque effect depend on the assumed mechanism 
for the initial energy deposition in the transverse plane 
and in space-time rapidity. Thus, the correlations in the 
longitudinal direction can serve as a tool to study the 
mechanisms of the initial energy deposition in the fire¬ 
ball. 

In our original study [7j, concerning the Au+Au col¬ 
lisions at the Brookhaven Relativistic Heavy-Ion Col¬ 
lider (RHIC), we have not carried out the event-by¬ 
event analysis, which is presented here. Instead, we 
have used the fact that the standard deviation of the 
torque angle between the forward and backward direc¬ 
tions in the Glauber initial condition has a particular 
value (of the order of 10° for the pseudorapidity separa¬ 
tion of 2 units). This value was implemented in the ini¬ 
tial shape, later evolved with hydrodynamics. We have 
also carried out the statistical hadronization at freeze- 
out with THERMINATOR [8! to generate the distribution 
of the observed hadrons. The finite number of the pro¬ 
duced hadrons leads to purely statistical fluctuations of 
the event-plane angles, which must be carefully sorted 
out through appropriate torque-effect measures. Ref. [7] 
also provided the estimates for two- and four-particle 
measures of the effect, as well as explored the central¬ 
ity dependence of the effect at the level of the Glauber 
initial state. 

In this work a full-fledged hydrodynamic event-by¬ 
event analysis of the torque effect is presented and com¬ 
pared to the results of the CMS Collaboration [91. Specif¬ 
ically, to prepare initial conditions we use the Glauber 
Monte Carlo model with the assumption of asymmet- 
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ric emission in the space-time rapidity from the target 
and projectile participants [TO], motivated with the pa- 
rameterizations from Refs. m- As mentioned above 
and discussed in detail in Ref. [7], the asymmetry of the 
energy deposition in pseudorapidity gives the torque of 
the event-plane as a function of the space-time rapidity. 
A similar mechanism yields the finite correlation range 
of particle multiplicities in pseudorapidity intervals as a 
function of the separation between the intervals m- We 
note that multiplicity correlation in different pseudora¬ 
pidity intervals have been the subject numerous experi¬ 
mental and theoretical investigations [12]. 

The event-plane decorrelation as a function of the 
pseudorapidity separation has been identified in several 
models of the initial conditions in relativistic heavy-ion 
collisions [13]. First experimental indications for the 
event-plane dependence on rapidity came from the az- 
imuthally sensitive Hanbury-Brown-Twiss analysis El- 
Very recently, the CMS Collaboration measured the ef¬ 
fect with the help of three-bin correlations [9j, which 
greatly reduces the sensitivity to non-flow effects. 


II. FORMALISM 

We perform event-by-event simulations for Au+Au col¬ 
lisions at snn = 200 GeV (RHIC) and Pb+Pb colli¬ 
sions at 2.76 TeV carried out at the Large Hadron Col¬ 
lider (LHC). To pinpoint the torque effect, we compare 
rapidity correlation for two sets of initial conditions: a 
forward-backward symmetric longitudinal profile of the 
deposited entropy for each wounded nucleon (no-torque, 
which serves as the reference calculation), and the profile 
with asymmetric deposition from the left- and right-going 
nucleons (torqued initial conditions). The subsequent ex¬ 
pansion of the fireball is modeled with 3+1-dimensional 
viscous hydrodynamics. Hydrodynamic simulations are 
carried out on event-by-event basis for every generated 
initial condition. The viscosity to entropy density ra¬ 
tio is ? 7 /s = 0.08 for the shear viscosity (independent of 
temperature) and £/s = 0.04 for the bulk viscosity (for 
T < 170 MeV) [15]. The fluid expands from the initial 
time Tg = 0.6 fm/c until freeze-out at the temperature 
Tp = 150 MeV. At freeze-out, hadrons are emitted statis¬ 
tically and afterwards resonance decays occur, as imple¬ 
mented in the event generator THERMINAT0R [8j. For each 
investigated centrality we generate from 200 to 300 hydro 
events. To increase statistics, from 200 to 2000 (depend¬ 
ing on the energy and collision centrality) THERMINAT0R 
events are generated for each freeze-out hypersurface. 

The initial entropy density in the fireball is calcu¬ 
lated with GLISSANDO pB, T7]. The positions Xi,yi of 
the N .|_ ( N _) wounded nucleons [T5] from the forward- 
(backward-)going nucleus are generated in each event. 
We use a Gaussian wounding profile for the nucleon- 
nucleon collisions in the Glauber model, and an excluded 
distance d = 0.9 fm for the nucleon configurations in 
the nuclei. The total inelastic nucleon-nucleon cross sec¬ 


tion is taken to be 42 mb at RHIC and 63 mb at the 
LHC. The entropy density in the transverse plane and 
the space-time rapidity Tyy is modeled as a sum of contri¬ 
butions from the individual nucleons 


JV+ AT_ 

s{x,y,V\\) = '^2gi(x,y)f+{r) + )+'^2g i (x,y)f_{ri ll ) . ( 1 ) 

i—1 i= 1 


The density in the transverse plane has a Gaussian form 
gi(x, y) = n [(1 - a) + A™ 11 a] x 

(x - Xi ) 2 + (y- yi) 2 


exp 


2a 2 


( 2 ) 


where the width of the Gaussian is a = 0.4 fm, N£° n is 
the number of collisions for the i-th nucleon, and a = 0.12 
and 0.15 and RHIC and the LHC, respectively. We note 
that E*= + i N ?° U = N i° n = N bin ~ the number of 

binary collisions. The non-zero value of a implements the 
mixed Glauber model El, adding binary collisions to the 
wounded nucleon model. The parameter k is adjusted for 
each energy to reproduce the average particle multiplicity 
in central events. 

The longitudinal density profile (an approximate “tri¬ 
angle” ) 


/±0?||) = Vb Z ±m m m ) for |?7||| < y heam (3) 

beam 


is a product of an asymmetric linear function and an over¬ 
all rapidity profile H{rj n), taken in the form of a plateau 
with Gaussian tails [2D], 


H{i 7||) = exp 


f (I 7 ?!! I — 7 7p) 2 Q(I 7 7|| I — Vp) 
\ 2u 2 


(4) 


where cr v = 1A and g p = 0.75 and 1.15 at RHIC and the 
LHC, respectively. The ansatz © with an asymmetric 
deposition of entropy is based on the analysis of asym¬ 
metric d-Au collisions at RHIC m- The event planes 
defined with the forward- and backward-going wounded 
nucleons lead to different event-plane directions, due to 
fluctuations of N+ and _/V_. The asymmetric deposition 
of entropy makes the event plane direction rapidity de¬ 
pendent [?'. 

We recall that the above formulas also lead to a tilt of 
the fireball away from the collisions axis (rotation in the 
z — x plane), which generates the rapidity-odd directed 
flow [ID]. The exact form of the asymmetric ansatz for 
the initial fireball can be varied, leading to slightly dif¬ 
ferent predictions for the directed flow mm- In the 
following we use the ansatz (J3|, which was shown to re¬ 
produce the directed flow at RHIC with the viscous hy¬ 
drodynamic evolution using g/s = 0.08 [22] , 

For comparison, we also perform simulations with a 
symmetric (no-torque, no-tilt) emission ansatz, where 
f±(? 7 ||) = H(i 7 ||). All the simulation points in this pa¬ 
per are drawn with error bars estimated from resampling 
of the Monte Carlo events obtained from THERMINAT0R. 
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FIG. 1. (color online) The 3-bin measure for the elliptic flow of charged hadrons for two centralities and for elliptic (n = 2) and 
triangular (n = 3) flow. The CMS data come from Ref. [9j. The empty squares correspond to the model without the torque 
effect and are compatible with 1. The filled squares are the result of simulations with the torque effect. The solid line shows 
the calculation with cumulative events (see text for details). The model points are slightly displaced horizontally for better 
visibility. 


III. 3-BIN CORRELATIONS 


We first show our results for the torque effect based 
on the three-bin correlation measure introduced by the 
CMS Collaboration j9j. It is defined as 


r n {v a ,ri b ) 


VnA(~V a 1 V b ) 

V n A (V a ,V b ) 


(5) 


where 


VnA(r] a , V b ) = ((cos [n((j> a - </> b )])). (6) 

The averaging is over all charged hadron pairs where one 
particle belongs to bin a and the other to bin b, and 
over events. The geometric interpretation of the above 
measure is very simple when the magnitude of the flow 
coefficient, v n (r / a ), is uncorrelated to the event-plane an¬ 
gle in bin a, 'F„(? 7 a ). In that case 


, a b , cos[n(T„,(-r ? a )-d' w (r ? fc ))] 
' ' ' cos[n(d'„(? 7 a ) — \E' n (ry b ))] 


(7) 


An important feature of measure ([5]) is that it cancels 
the random fluctuations of the event-plane angles, which 
result from a finite number of hadrons in the pseudo¬ 
rapidity bins. For distant bins a and 6, the fluctuation 
of (’F„(? 7 a ) — d'„(r/ b ) is on the average the same as for 


(H/ n (— rj a ) — 'f’n('n b )- Then, in the absence of a genuine 
torque effect, we would have r n (rj a , r] b ) = 1. The de¬ 
parture of the measure from unity quantifies the torque 
effect. 

Our numerical simulations are carried out for the same 
bin definition and kinematic cuts similar to Ref. [5]: 
-2.4 < rf < 2.4, 4 < r] b < 5, 0.3 < p% < 3 GeV, 
0 < p b T < 3 GeV. The results are shown in Fig. [l] for 
two sample centralities (0-5% and 30-40%) and for n = 2 
(elliptic flow) and n = 3 (triangular flow). First, we note 
that while the statistical errors are large due to limited 
statistics, one may clearly see a qualitative agreement 
of the model points with the experiment. Also, the no¬ 
torque simulation (empty squares) is compatible with 1 
within the statistical uncertainties, in agreement with the 
discussion above. 

To circumvent the limits of the finite statistics, we have 
also carried out the analysis with cumulative events (solid 
lines in Fig. [l]). In this case we evaluate by summing 
up hadrons from all THERMINATOR events simulated on 
top of a given hydro event. That way we eliminate the 
fluctuations coming from a finite number of the produced 
particles and minimize the non-flow effects. This trick is 
possible in model calculations, as opposed to the exper¬ 
iment, since in each THERMINATOR event run on top of 
a given hydrodynamic event we have a fixed geometry 
of the collision. We note that the result of the analysis 
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FIG. 2. (color online) Model predictions for the torque measure rij -13111 for charged particles at the LHC (panels a and c) and 
RHIC (panels b and d) energies and for two sample centralities. The empty squares correspond to the model without the 
torque effect. The filled squares are the result of simulations with the torque effect. The solid line shows the calculation with 
cumulative events. 


with the cumulative events yields a smaller torque than 
the experiment. 


IV. 2-BIN CORRELATIONS 


Next, we pass to the results for the torque measure 
introduced in Ref. [71. It is defined via two-particle cu- 
mulants with particles belonging to bins F and B: 


r 


2—bin 
2 




( 8 ) 


This measure is sensitive to the non-flow effects as well, 
hence is useful to estimate such contributions to the 
torque effect. 

The results of our simulations are shown in Figs. [2] and 
[3] for the cases of the elliptic and triangular flow, respec¬ 
tively. The cuts are 0.15 < pr < 3 GeV for RHIC and 
0.3 < pt < 3 GeV for the LHC. The corresponding 
calculation with cumulative events is indicated with the 
solid line. We note that due to non-flow effects the re¬ 
sults for the case without the torque (empty squares) are 
significantly below unity. The presence of torque (filled 
squares) lowers further the results, but this effect would 
be difficult to assess experimentally. We note that the 
difference of the torque and no-torque curves is very close 
to the corresponding difference for the cumulative events 
(where the no-torque result overlaps with unity). We 


also note that the predicted effect of the torque is signif¬ 
icantly higher at RHIC than at the LHC, and is stronger 
for n = 3 than for n = 2. 

Comparison of the results of Fig. [I] with Figs. [2] and [3] 
shows that the 3-bin measure is advantageous with regard 
to the 2-bin measures, as the non-flow effects are greatly 
reduced. 


V. CONCLUSIONS 

The very accurate studies at the LHC make the flow 
correlation studies feasible and practical. Indeed, the er¬ 
ror bars in the CMS analysis of Ref. [9j are much smaller 
than the size of the investigated effect. Our event-by¬ 
event simulations confirm that the interpretation of the 
behavior of the CMS measure r n (r] a ,rf > ) via the torque 
effect proposed in Ref. [7] is natural: the decorrelation of 
the event-plane orientations between bins distant in pseu¬ 
dorapidity follows from fluctuations of the left- and right¬ 
going wounded nucleons and their asymmetric emission 
profiles. 

Nevertheless, precise modeling of the torque effect may 
be a numerical challenge, as it requires a very large event- 
by-event statistics. Yet, with the number of events at our 
disposal, we can clearly see several features: the effect is 
larger at RHIC than at the LHC, it is also larger for the 
triangular flow than for elliptic flow. The presence of non- 
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FIG. 3. (color online) Same as Fig. [ 2 ] but 


flow affects the two-bin measures, while the application of 
the three-bin method proposed by the CMS J9] minimizes 
these effects. 

In future studies one could also apply to our hydrody¬ 
namic model the Principal Component Analysis (PCA), 
recently brought up in Refs. [231 M- This would allow 
one to quantify the torque effect (and other correlation 
effects) from a more general perspective of the multivari¬ 
ate analysis, similar to the case of the AMPT model [23] 
investigated in Ref. [23] [24] • In analogy to the two-bin 
measure r 2 _bm , dominated by non-flow effects, we expect 
that these effects will also influence the PCA technique 
based on covariance matrix composed of two-bin observ¬ 
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the triangular flow torque measure -bm 


ables. Thus it is important to compare model predictions 
with and without the torque effect in the fireball fluctu¬ 
ations, as done in the present work. 
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